A new fuzzy rule based multi-objective optimization method for cross-scale injection molding of protein electrophoresis microfluidic chips

Injection molding is one of the most promising technologies for the large-scale production and application of polymeric microfluidic chips. The multi-objective optimization of injection molding process for substrate and cover plate on protein electrophoresis microfluidic chip is performed to solve the problem that the forming precision is difficult to coordinate because of the cross-scale structure characteristics for chip in this paper. The innovation for this research is that an optimization approach and a detailed fuzzy rule determination method are proposed in multi-objective optimization for protein electrophoresis microfluidic chip. In more detail, firstly, according to the number and level of process parameters, the orthogonal experimental design is carried out. Then, the experiments are performed. Secondly, the grey relational analysis (GRA) approach is employed to process the response data to gain the grey relational coefficient (GRC). Thirdly, the grey fuzzy decision making method which combines triangular membership function and gaussian membership function is adopted to obtain the grey fuzzy grade (GFG). After that, the optimal scheme of process parameters was predicted by the grey fuzzy grade analysis. Finally, the superiority of Taguchi grey fuzzy decision making method are verified by comparing the results of original scheme, optimal scheme and prediction scheme. As a result, compared with the original design, the residual stress of substrate plate (RSS), residual stress of cover plate (RSC), warpage of substrate plate (WS), warpage of cover plate (WC) and replication fidelity of microchannel for substrate plate (RFM) on the prediction scheme for Taguchi grey fuzzy decision making method were reduced by 32.816%, 29.977%, 88.571%, 74.390% and 46.453%, respectively.

. Due to the narrow process window and many parameters, the problem of how to obtain precise chip forming process efficiently has not been solved yet. During the forming process, the cross-scale geometric deviation between the macroscopic warping deformation of chip and the deformation of microchannel is difficult to coordinate, and the accuracy is difficult to ensure. Meanwhile, a qualified microchannel chip needs not only high replication fidelity of microchannel, but also good bonding quality. However, the warpage and residual stress of substrate and cover plate for microfluidic ship will seriously affect the bonding quality of the chip, which will lead to non-sealing and deformation of the chip. Most of the existing researches only consider one objective to carry out single-objective optimization [16][17][18] . Nevertheless, the residual stress, the warpage and the replication fidelity of microchannel for the substrate and cover on microfluidic chips should be considered simultaneously to meet the needs of high-performance chips. Few studies have explicitly proposed a detailed method for multi-objective optimization of injection molding on microfluidic chip. Similar to the optimization of all Microinjection molding products, there are two primary assessment criterions on microfluidic chip molding: the shape accuracy of plate and the molding quality of microchannel, which are in conflict to some extent (the plates of substrate and the cover being formed in the same mold simultaneously). Therefore, it is necessary to establish a balance between them. The optimization of substrate and cover plates for microfluidic chip needs to consider the influence of multiple aspects on the overall performance (the replication fidelity of microchannel, the warpage and the residual stress, etc.), which is a multi-objective optimization problem (MOOP). In the past three decades, scholars have put forward many effective approaches to solve this problem. Normally, there are two main directions employed to deal with the multi-objective problem. One is the global optimization method combining optimization algorithm and surrogate model, named SMOA; the other is the method using a multi-criterial decision method to search the best solution among all sample points, named MCDM (Table 1). To be more specific: (1) In order to obtain the global optimal solution, many optimization algorithms are proposed, such as: genetic algorithm (GA) 19 , non-dominated sorting genetic algorithm (NSGA-II) 20 , multi-objective particle swarm optimization algorithm (MOPSO) 21 , multi-objective genetic algorithm (MOGA) 22 , etc. (2) In order to reduce the experimental workload and improve the optimization efficiency, a large number of surrogate models are adopted. Such as: Kriging model 23 , radial basis function (RBF) 24 , response surface methodology (RSM) 25 , etc. (3) To quickly determine the best compromise, multi-criterial decision making (MCDM) methods are also widely employed. Such as: grey relational analysis (GRA) 26 , technique for order preference by similarity to ideal solution (TOPSIS) 27 , fuzzy logics (FL) 28 etc.
According to the literature on multi-objective optimization of injection molding, the multi-objective optimization based on surrogate model and optimization algorithm is the mainstream optimization method at present. However, most of the existing algorithms can work effectively with less than three optimization objectives. In SMOA, when more optimization objectives need to be considered, the application of optimization algorithms will be limited 29,30 . In addition, in order to obtain an accurate surrogate model, a large number of sample points are needed, which is a great challenge to establish an accurate surrogate model. Moreover, because of the strong nonlinearity of the replication fidelity for microchannel on microfluidic chip, the low precision of the established surrogate model will lead to a large deviation between the results of the optimal scheme obtained by the optimization algorithm and the experimental value, which can not ensure the effectiveness of the global optimal scheme (as shown in Table 7). Therefore, surrogate model and optimization algorithm is not suitable for multi-objective optimization in injection molding of microfluidic chip.
In MCDM, the GRA is also widely applied to solve multi-objective optimization problems. Therefore, a method combining GRA and Taguchi design was proposed by researchers 31,32 , which can not only overcome the defect that Taguchi method cannot solve multi-objective optimization 33 , but also mine the information of the whole parameter space with limited sample size as much as possible. However, the GRA cannot provide an optimal solution of highly robust for a given multi-objective optimization problem, because this method cannot quantitatively or qualitatively distinguish the ideal case of foraging problem with no solution (black) and unique solution (white) 34 . In detail, different normalized formulas are applied on the basis of target characteristics, and then the grey relational coefficient and GRG are calculated, which leads to the uncertainty of the optimal solution of GRA 35,36 . Thus, applying GRA method to solve multi-objective optimization problems has certain limitations.
Fuzzy logic was proposed by Zadeh in 1965, which is usually adopted to solve problems of uncertainty and ambiguity 37 . In addition, fuzzy logic has been applied to the field of multi-objective optimization in recent  40 . At present, the research of multi-objective fuzzy decision making is relatively simple, and there is no definition and guidance for the writing of fuzzy rules. In addition, the existing multi-objective fuzzy decision-making system only applies to the case where the weights of each index are equal, and there is no effective scheme for the case where the weights of each index are different. Therefore, a novel fuzzy decision making method and specific fuzzy rule writing method are proposed to solve the problem of different objective weights and rule generate in multi-objective optimization. In this paper, the Taguchi design, the grey relational analysis and fuzzy decision method are simultaneously applied to optimize the injection molding process parameters of substrate and cover plates for typical microfluidic chip (protein electrophoresis microfluidic chip).

Protein electrophoresis microfluidic chip experiment
Design for protein electrophoresis microfluidic chip. In this paper, a cross protein electrophoresis microfluidic chip with a single channel structure was designed. The chip is composed of substrate plate and cover plate, of which the length, width and thickness for the substrate plate are 50 mm, 28 mm and 0.8 mm respectively, while the length, width and thickness for the cover plate are 50 mm, 28 mm and 0.6 mm respectively. The detailed structure and size of chip are shown in Fig. 1.
Experimental material and equipment. PMMA was selected as the material for the production of protein electrophoresis microfluidic chips due to its advantages of high light transmittance, good solvent and chemical compatibility in this study. To be more specific, the PMMA CM-205 produced by ChiMei from Taiwan in China was applied, and the properties are shown in Table 2. In addition, the injection molding system (Fig. 2a) includes three parts: the precision injection molding machine, the feeder and the mold temperature machine. This molding system has the advantages of programmable, stable process performance and high temperature control accuracy. The detailed product models, company and production areas are shown in Table 3. Meanwhile, a dual-type cavity mold core that can simultaneously form the substrate and cover plate for the chip was adopted, as shown in Fig. 2b.

The theory of multi-objective optimization for protein electrophoresis microfluidic chip
The multi-objective optimization method. In this paper, the Taguchi orthogonal design method, the grey relational analysis method, the grey Taguchi fuzzy decision method and the factor influence analysis are simultaneously adopted in multi-objective optimization of protein electrophoresis microfluidic chip. Firstly, the design approach for multi-objective optimization is determined and the optimization objectives are determined according to the application requirements of protein electrophoresis microfluidic chip. Secondly, the change level of design variables is determined according to the original design parameters. Then, the orthogonal experimental design is carried out. Thirdly, the response data of optimization objectives is obtained by performing  www.nature.com/scientificreports/ practical experiments. Fourthly, the grey relational analysis method is applied to process the response data to obtain the grey relation grade. In addition, the membership functions for input variables, the membership functions for output variables and corresponding fuzzy subset levels are confirmed. Meanwhile, the fuzzy rules are compiled according to the weights of optimization target. Moreover, the grey fuzzy grade (GFG) was analyzed to predict the optimal combination of process parameters. Finally, the superiority and effectiveness of Taguchi grey fuzzy decision making method are verified by comparing the original scheme, the optimal scheme and the prediction scheme. The flow chart of multi-objective optimization method for protein electrophoresis microfluidic chip is shown in Fig. 3.
The evaluation criteria of multi-objective optimization. The optimization of protein electrophoresis microfluidic chip has two main evaluation criteria: the shape accuracy of plate and the molding quality of microchannel. To be more specific, the protein electrophoresis microfluidic chip will warp and crack during the storage process because of the residual stress. Hence, the residual stress is one of the shape accuracy criteria. Meanwhile, warpage directly affects the bonding quality of protein electrophoresis microfluidic chips. Therefore, the warpage is one of the shape accuracy criteria. As the main functional unit of protein electrophoresis microfluidic chip, the molding quality of microchannel is particularly important. Accordingly, the replication fidelity is the molding quality of microchannel for evaluation criteria. In particular, the residual stress of the protein electrophoresis microfluidic chip can be characterized by the birefringence effect of the chip. The relationship between refraction value and residual stress can be expressed as σ = �S K×H , where σ is the residual stress, S represents the refraction value, K donates the pressure optical coefficient, H is the thickness of the measured part. The larger the residual stress, the stronger the birefringence effect and the larger the measured value. Likewise, the warpage of the protein electrophoresis microfluidic chip can be characterized by the distance between the maximum value and the minimum value in the vertical direction of the chip surface. In addition, the replication fidelity of microchannel is represented by the mean square root of the ideal microchannel contour and the actual microchannel contour. The formula is as follows 41 : where MSR represents the value of the mean square root, n i is the point of the actual microchannel contour, N i denotes the point of the ideal microchannel contour, i is the amount of the point for microchannel contour. Table 4 and Fig. 4 show the test equipment and models adopted for residual stress, warpage and replication fidelity of microchannel for protein electrophoresis microfluidic chip.
The evaluation criteria for protein electrophoresis microfluidic chip include: (1) the value of residual stress for substrate plate and cover plate, (2) the value of warpage for substrate plate and cover plate, (3) the replication fidelity of microchannel for substrate plate. In this study, the birefringence value for substrate plate and cover plate represents the residual stress for substrate plate and cover plate (Fig. 5a). In Fig. 5a, the substrate plate and cover plate of the microfluidic chip manufactured by injection molding are shown on the left of figure, and the Figure 2. The injection molding system and the mold core of protein electrophoresis microfluidic chips (a) The injection molding system (b) The mold core. Table 3. The product models and production areas of injection molding system. www.nature.com/scientificreports/ birefringence effect of the substrate plate and cover plate is shown on the right of figure. The distance between the maximum and minimum value of line B in transverse symmetry for substrate and cover plate represents the warpage for substrate plate and cover plate (Fig. 5b). The mean square root of the microchannel at point A, 20 mm from the short edge of the substrate, denotes the replication fidelity of microchannel for substrate plate (Fig. 5c). In Fig. 5c, the positions of microchannels to be measured in the substrate are shown on the left side of the figure. The microchannel morphology of point A measured by laser confocal microscopy is shown in the middle of the figure. The comparison between the actual and ideal microchannel morphology is shown in the right side of the figure.

The design of experiment design for multi-objective optimization
Design variables. According to the existing researches on the replication fidelity of microchannel, the warpage and the residual stress of microfluidic chips in injection molding 14,[16][17][18] and experience in protein electrophoresis microfluidic chip production, five major process parameters were selected as design variables, including the melt temperature (MT), the injection pressure (IP), the injection speed (IS), the packing pressure (PP) and the packing time (PT), in this paper. Five experimental levels are selected for all design variables. The range of process parameters on injection molding for protein electrophoresis microfluidic chip is shown in Table 5.
The design and results of experiment. An Orthogonal experimental design, L25(5 5 ), was carried out for the injection molding process parameters of five factors and five levels in Table 5. Then, according to orthogonal experimental design, the practical experiments are carried out, as shown in Table 6 (each experiment result is the average value of three experiments under the condition of stable process).   www.nature.com/scientificreports/ Optimization objectives. In the process of protein electrophoresis microfluidic chip optimization, various criteria of substrate plate and cover plate should be considered. In general, the residual stress and the warpage are employed to assess the shape accuracy. In addition, the replication fidelity of microchannel is applied to evaluate the molding quality of microchannel. Hence, five optimization objectives are selected for multiobjective optimization of protein electrophoresis microfluidic chip in this research, which are the residual stress of substrate, the residual stress of cover, the warpage of substrate, the warpage of cover and the replication fidelity of microchannel. The specific objectives are as follows: (1) The value of birefringence for substrate plate; (2) The value of birefringence for cover plate; (3) The distance between the maximum and minimum value of line B in transverse symmetry for substrate plate; (4) The distance between the maximum and minimum value of line B in transverse symmetry for cover plate; (5) The mean square root of the ideal microchannel contour and the actual microchannel contour for the microchannel at point A on the substrate.

Taguchi grey fuzzy decision making method
Grey relational analysis method. According to the experimental results of protein electrophoresis microfluidic chip shown in Table 6, the process parameters of design variables MT, IP, IS, PP and PT are taken as the input of the model, while the RSS, RSC, WS, WC and RFM are the output. Meanwhile, different methods were adopted to fit the surrogate models [Kriging model, FRB model, RSM (first order), RSM (second order) and RSM (third order)]. The accuracy (R 2 ) of each surrogate model is shown in Table 7. As can be seen from Table 7, the highest accuracy is 0.880 and the lowest accuracy is 0 in the surrogate model. Because of the low  Grey relational analysis is a multi-index decision making method proposed by Deng in the 1980s 34 . The GRA can transform multi-objective optimization problems into single-objective optimization problems to solve the disadvantages and limitations of Taguchi method that can only deal with single-objective problems. In the data   www.nature.com/scientificreports/ preprocessing stage, different normalization formulas are selected to deal with the original sequences according to the different directions of optimization objective characteristic. To be more specific, on the basis of the characteristics of the target direction "the-bigger-the-better", "the-lower-the-better", and "the-closer to the target objective value-the-better", three normalization formulas can be adopted in this approach.
If the optimized direction of original sequence is "the-lower-the-better", the original sequence can be normalized as follows 30 : If the optimized direction for original sequence is "the-bigger-the-better", the original sequence can be normalized as follows: If the optimized direction of original sequence is a target value, namely the feature for original sequence is "the-closer to the target objective value-the-better", the original sequence can be normalized as follows: where x * i (t) represents the sequence generated by the grey relation, x i (t) denotes the original sequence of experiment results for ith element in the sequence of tth, Max t x i (t) is the maximum value for the sequence of tth, Min t x i (t) represents the minimum value for the sequence of tth, i = 1, 2, 3…n 1 and t = 1, 2, 3,…n 2 , n 1 is the sample size of Taguchi design and n 2 denotes the quantity of optimization objectives, T is the target value.
In this paper, the residual stress of substrate plate, the residual stress of cover plate, the warpage of substrate plate, the warpage of cover plate and the replication fidelity of microchannel should be lowest. Hence, formula (2) is applied to normalize the original sequences of experiment results to [0, 1], as shown in Table 8. Normally, the reference sequence is defined as 1, representing the optimum performances in theory. Therefore, the comparable sequence which is closest to 1 will be considered the optimal scheme.
The relation between reference sequence and comparable sequence can be confirmed by grey relational coefficient (GRC). Then, the calculation formula is as follows: where γ x * r (t), x * i (t) is the grey relational coefficient, u denotes the distinguishing coefficient, u ∈ [0,1], in general, the value of distinguishing coefficient is 0.5, � ri (t) , min , max can be deduced by: Formula (5) is adopted to calculate the normalized results to obtain GRC ( Table 8). The grey relational grade (GRG) can be gained by employing GRC. If the weight for each optimization objective is equal, the following formula is used for calculation: If the weight of each optimization objective is different, GRG can be calculated as follows: .
. www.nature.com/scientificreports/ where w t denotes the weight of tth optimization objective, n represents the amount of the optimization objectives. In this study, the weights of optimization objectives are different. Hence, GRG is obtained by using formula (10). Meanwhile, the importance of each optimization objective is the replication fidelity of microchannel, the warpage and the residual stress. Moreover, the weights for warpage on the plates of substrate and cover are the same, as is the residual stress for the plates of substrate and cover. The weights of the residual stress of substrate plate, the residual stress of cover plate, the warpage of substrate plate, the warpage of cover plate and the replication fidelity of microchannel are set as 0.1250, 0.1250, 0.1750, 01,750 and 0.4000, respectively ( Table 9). As shown in Table 10, the GRGs of experimental results are ranked from big to small, while the No.3 shows the biggest value of GRG in 25 experiments. Therefore, the No.3 is the best scheme for multi-objective optimization of protein electrophoresis microfluidic chip.
Grey Fuzzy decision making method. Fuzzy logic, first proposed by Zadeh, is a method to deal with imprecise and undefined boundary problems 36 . Fuzzy decision systems can reason output and calculation results   www.nature.com/scientificreports/ by using fuzzy set theory, and fuzzy sets can provide accurate description by using mathematical expressions. in this study, in order to solve the problems of the uncertainty in GRA optimization results, the fuzzy decision method under different weights and the writing of fuzzy rules, a new fuzzy decision method and specific writing of fuzzy rules are proposed. The new grey fuzzy decision making method can be divided into four steps for analysis and calculation. Meanwhile, the software of MATLAB ® is adopted for fuzzy logic analysis in this study.
Step 1 Fuzzification. The principle of fuzzification is to apply language variables to transform clear values into fuzzy quantities. Different and appropriate membership functions are selected to assign membership degrees for each language item. There are many membership functions that can be applied, for instance: triangular function, gaussian function, trapezoidal function, γ function and sigmoidal membership function, etc 38 . Among these membership functions, the triangular function was widely adopted due to its simplicity and high calculation efficiency in the existing study 42 . However, in the case of many optimization objectives and high uncertainty of optimization scheme, the selection of triangular function will cause different schemes to have the same fuzzy grade and the order of optimization results is not clear (Table 16). Therefore, the triangular function and the gaussian function are adopted simultaneously to avoid the situation of the same fuzzy grade for output in this study, namely, the gaussian function is employed for input and the triangular function is employed for output. The triangular and gaussian functions are defined as follows 34,37,38 : where x represents the variable; a, b, and c are the vertices of a triangle; σ denotes the coefficient of a gaussian function; m is the center of symmetry of a gaussian function. In this paper, the GRCs of the optimization objectives are adopted as the input variables. Then, the gaussian membership function (12) is applied to fuzzify GRCs into 3 levels: high (H), medium (M) and low (L), namely, the fuzzy subsets of input are divided as 3 grades: H, M and L, as shown in Fig. 6. Correspondingly, the output variables are defined as 11 levels using the triangular membership function (11), which are very very very high (VVVH), very very high (VVH), very high (VH), high (H), between medium and high (MH), medium (M), between medium and low (ML), low(L), very low (VL), very very low (VVL) and very very very low (VVVL) respectively. Likely, the fuzzy subsets of output are divided as 11 grades: VVVL, VVL, VL, H, MH, M, ML, L, VL, VVL and VVVL, as shown in Fig. 7. The value range of fuzzy subset of input and output is shown in Table 11.
Step 2 Fuzzy rule base. Fuzzy rule base is a collection of fuzzy rules, which can be applied to express the relationship between input variables and output variables. In fuzzy decision making, fuzzy rules are the key to perform reasoning. At present, IF-THEN rule is most commonly used and can be expressed as: The new and detailed grey fuzzy rules are written as follows: ① Carried out the full-factor design according to the optimization objectives and fuzzy subset grades of the optimization objectives (the fuzzy subset grades of input). ② Divide and match the value 1 according to the grade of the output fuzzy subsets. ③ Redistribute if premise (antecedent), then conclusion consequent . To be more specific, during this research, the fuzzy rule base is defined by the criterion that the larger the GRC of each optimization objective is, the higher the GFG value is. the input of fuzzy rule base is composed of five optimization objectives, which are RSS, RSC, WS, WC and RFM respectively. Meanwhile, one output is set, which is GFG. The number of fuzzy subset grades for optimization objectives is 3, and the number of optimization objectives is 5. Therefore, the amount of fuzzy rules should be 243 (3 5 ). The full-factor design based on optimization objectives (RSS, RSC, WS, WC and RFM) and fuzzy subset grades of input (H, M and L) is shown in Table 12. The weight values of the fuzzy subset grades on the input for each optimization objective are shown in Table 13. Namely, the weights of H, M and L grades on the fuzzy subset for RSS and RSC are 0.1250, 0.0625 and 0, respectively; the weights of H, M and L grades on the fuzzy subset for WS and WC are 0.1750, 0.0875 and 0, respectively; the weights of H, M and L grades on the fuzzy subset for RFM are 0.4000, 0.2000 and 0, respectively. The values and the fuzzy subset grades of output are shown in Table 14. Namely, the corresponding values the VVVH, VVH, VH, H, MH, M, ML, L, VL, VVL and VVVL grade for the output fuzzy subsets are 1.0, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 and 0, respectively. The weight values of the grades for each input fuzzy subset, the comprehensive weight value and the grades of corresponding output fuzzy subsets are shown in Table 15 (the "CV" column is the comprehensive weight). The final fuzzy rule scheme of this study is as follows: Rule 1 : if "RSS is H and RSC is H and WS is H and WC is H and RFM is H"then "GF is VVVH"; else  www.nature.com/scientificreports/ Rule 2 : if "RSS is L and RSC is L and WS is H and WC is M and RFM is H"then"GFGisH"; else · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · (13) Rule 243 : if "RSS is L and RSC is L and WS is L and WC is L and RFM is L"then"GFGisVVVL" Table 12. The full-factor design based on optimization objectives and fuzzy subset grades of input.

No RSS RSC WS WC RFM RSS RSC WS WC RFM
Step 3 Fuzzy inference. The maximum and minimum operation on Mamdani approach is adopted to execute fuzzy inference for multiple response outputs, and the fuzzy output values are obtained by using fuzzy rule base. The output formula of fuzzy inference is: where the F i is the fuzzy subsets of output defined by membership functions; y i denotes the value of output; x i represents the value of input for mth optimization objective; β is the range of coefficient; A Xm denotes the mth fuzzy rules; βF i y i is the grade of fuzzy inference for output. In this paper, formula (14) is employed to calculate the grade of grey fuzzy inference.
Step 4 Defuzzification. The principle of defuzzification is to convert the language output into clear and concrete values. Currently, the center of gravity method can effectively transform fuzzy inference into clear value, which is called fuzzy grade as y i . The calculation formula is as follows: In this study, (15) is adopted to convert the calculated grades of grey fuzzy inference into clarity values.

Results and discussion
Grey fuzzy grade. The GFGs combining the triangular membership function and gaussian membership function, the GFGs of the triangular membership function and the GRGs of the GRA method are listed and arranged from large to small in Table 16. As can be seen from Table 16, some GFGs obtained by adopting triangular membership function are the same, which leads to unclear ordering of optimization results (No. 1, 3, 5 and 8 have the same GFGs). However, this problem does not exist in the combination of triangular membership function and gaussian membership function. Therefore, the GFGs obtained by combining triangular membership function with gaussian membership function has higher accuracy and discrimination. By comparing the GFGs combining the triangular membership function and gaussian membership function with the GRGs of the GRA method, it can be found that the rank of GFGs and GRGs are mostly the same, except for a few schemes with some differences. Meanwhile, compared with the same scheme, the value of GFG is greater than the GRG. Figure 8 makes a more intuitive comparison between the GFGs combining the triangular membership function and the gaussian membership function and the GRGs of the GRA method. The uncertainty of the scheme with grey fuzzy decision method is obviously reduced, and the grade is closer to the ideal value 1. Therefore, the method of Taguchi grey fuzzy decision making can provide a more reliable, qualitative and robust multi-objective optimization platform. In addition, as shown in Table 16, the value of GFG (the triangular membership function and gaussian membership function) and GRG (GRA) of scheme No.3 is the largest. Thus, the No.3 is the optimal compromise among 25 experiments.
Influence analysis of factors. In Taguchi grey-fuzzy decision making, the average grey fuzzy grade (AGFG) on each design variable for the same level is a significant target to decide the importance on design variables and the optimal parameter combination. The AGFG on each design variable for the same level can be obtained by classification and calculation of each design variable in orthogonal table at the same level. The calculation formula is as follows: where the v ou represents the AGFG of oth design variable on uth level; N ou is the GFG of oth design variable on uth level; o denotes the number of levels; u is the number of design variables.
In this study, formula (16) is employed to classify and calculate the AGFGs (the triangular membership function and the gaussian membership function) in Table 16. The calculation results are shown in Table 17. Moreover, the main effect diagram (Fig. 9) clearly shows the influence of design variables on each optimization objective. It can be seen from Table 17 and Fig. 9 that the maximum values of AGFGs for design variables MT, IP, IS, PP and PT are 0.8182, 0.7458, 0.7378, 0.7692 and 0.7462, respectively. Therefore, the scheme of MT 1 , IP 2 , IS 4 , PP 3 and PT 2 (MT: 250 °C, IP:160 MPa, IS: 65 cm 3 /s, PP: 80 MPa and PT: 5 s) is the optimal parameter combination that can simultaneously improve the performance of each optimization objective.
As shown in Table 17, the differences value between the minimum and maximum of AGFGs, arranged from big to small, are the TM, PP, IS, IP and PT. Hence, the design variable (MT) is the most significant to improve the performance of each optimization target for protein electrophoresis microfluidic chip. In addition, all design variables have uncertain effects on the performance of each optimization objective for protein electrophoresis microfluidic chip, according to the slope of GFG curve in Fig. 9. This results also show that the multi-objective optimization of process parameters is necessary to ensure the integrated performance for protein electrophoresis microfluidic chips.
Optimal scheme prediction and validation. In order to verify the feasibility and performance improvement effect of the optimal scheme for protein electrophoresis microfluidic chip, the GFG of the optimal scheme must be predicted theoretically and tested in practice after determining the optimal scheme of protein electrophoresis microfluidic chip. The prediction formula of GFG is 38 : www.nature.com/scientificreports/ where τ p is the total mean value of the AGFGs; τ h denotes the mean value of the best parameter combination in the AGFGs; u is the amount for design variables. The experimental results and the GFG of the prediction scheme, the optimal scheme and the original scheme are shown in Table 18 (each experiment result is the average value of three experiments under the condition of stable process). Specifically, the "Original" column is the original design, the "Triangular + gaussian" column is the grey fuzzy decision making method applying the triangular membership function and the gaussian membership function, the "Optimal" column is the best scheme obtained by adopting the grey fuzzy decision making method with the triangular membership function and the gaussian membership function, the "Prediction" column is the method of prediction, the "Theory" column is the theory value of prediction method, the "Experiment" column is the experiment value of prediction method, the "Error (%)" column is the relative deviation between the prediction result and the optimal result, the "Improvement (%)" column is the relative deviation between the prediction result and the original design, the "Grey fuzzy grade (new)" column is the grey fuzzy grade obtained by grey fuzzy decision making method after adding the prediction scheme and original design (27 sample points), the "Grey fuzzy grade (old)" column is the GFG obtained by grey fuzzy decision making method (25 sample points). It can be seen from Table 18 that the GFG of the theory prediction (1.008) and the GFG of the experiment prediction (0.902) are differing greatly. The reason for this difference is that after the addition of the prediction scheme and the original scheme, the minimum value of the target sequence is changed, which leads to the change of all the grey relational coefficients and the change of the overall GFGs. (the difference between the Grey fuzzy grade (new) and the grey fuzzy grade (old)). In addition, compared with the original design, the RSS, RSC, WS, WC and RFM of the prediction scheme were reduced by 32.816%, 29.977%, 88.571%, 74.390% and 46.453%, respectively. Compared with the optimal scheme of grey fuzzy decision making method, the RSS, RSC, WS, WC and RFM of the prediction scheme were reduced by − 12.016%, − 5.192%, 12.195%, 61.111% and 30.210%, respectively. Moreover, the GFG of the original design, optimal scheme and prediction scheme is 0.616, 0.874 and 0.902, respectively. The GFGs ranking (high to low) are the prediction scheme, optimal scheme and original design. Figure 10 shows a comparison result of the original scheme, the optimal scheme and the prediction scheme in an experiment. In conclusion, the optimal scheme is the prediction scheme using grey fuzzy decision making method. Accordingly, Taguchi grey fuzzy decision making method can be adopted to optimize protein electrophoresis microfluidic chip effectively.  www.nature.com/scientificreports/

Conclusion
During this research, a multi-objective optimization approach and detailed optimization process are proposed for substrate and cover production on protein electrophoresis microfluidic chip, and the effectiveness of the prediction scheme is evaluated through experiments. In more detail, the Taguchi orthogonal design method, the grey relational analysis method, the Taguchi grey fuzzy decision making method and the factor influence analysis are simultaneously adopted in multi-objective optimization of protein electrophoresis microfluidic chip.
The main conclusions of this study include: (1) In the optimization of protein electrophoresis microfluidic chip, the grey fuzzy decision making method using the triangular membership function has the problem of insufficient accuracy, which will lead to the same GFG in different schemes, while the grey fuzzy decision method which combines triangular membership function and gaussian membership function can solve this problem and get the appropriate GFG. (2) The detailed fuzzy rule determination method can effectively solve the multi-objective optimization problem with different weights of optimization objectives. (3) Compared with the GRG of the GRA method, the GFG of the grey fuzzy decision making method is closer to the ideal value. The grey fuzzy decision making method can provide a more reliable, qualitative and robust multi-objective optimization platform. (4) The results show that the optimal scheme is achieved when the MT is 250 °C, the IP is 160 MPa, the IS is 65 cm 3 /s, the PP is 80 MPa and the PT is 5 s. (5) Compared with the original design, the RSS, RSC, WS, WC and RFM of the prediction scheme for Taguchi grey fuzzy decision making method were reduced by 32.816%, 29.977%, 88.571%, 74.390% and 46.453%, respectively.

Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].